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Abstract 

The modest aim of this short article is to provide some new results for a screw 
dislocation in a functionally graded material within the theory of gradient elastic- 
ity. These results, based on a displacement formulation and the Fourier transform 
technique, extends earlier findings (Lazar, M., 2007. On a screw dislocation in a 
functionally graded material, Mech. Res. Comm. 34, 305-311) obtained by the 
stress function method, to the case of the second gradient elasticity. Rigorous and 
easy-to-use analytical expressions for the displacements, the strains and the stresses 
are obtained which are free from singularities at the dislocation line. 

keywords: Screw dislocation; Functionally graded material; Second strain gra- 
dient elasticity 

1 Introduction 

The stress function technique was employed by Lazar [1] to derive stress and strain 
fields for a screw dislocation in a functionally graded material by using the theory 
of first gradient elasticity. Analytical non-singular expressions were derived for the 
strains and the (first order) stresses, but the double stresses remained singular. Such 
results are derived here by using the Fourier transform technique which, in addi- 
tion, provides exact analytical expressions for the displacement field. Moreover, the 
problem is reconsidered within the framework of "second strain gradient elasticity" 
which eliminates the singularities from the double stress expressions as well. Recent 
work on gradient elasticity [2,3] has revealed the need of using higher order gradi- 
ents of strain in the stress-strain relation in order ro interpret experimental results 
pertaining to dislocation density tensor and more accurately describe the details 
of the relevant stress/strain fields near the core of dislocations contained in small 
volumes. This is the case in particular, for dislocations contained in functionally 
graded materials (FGMs), the use of which has advocated since mid 80's in relation 
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to ultra high temperature and ultra high weight requirements for aircraft, space 
vehicles and other applications. Generally FGMs refer to heterogeneous composite 
materials, in which mechanical properties are intentionally made to vary smoothly 
and continuously from point to point. This is controlled by the variation of the vol- 
ume fraction of the constituent materials. Ceramic/ceramic and metal/ceramic are 
typical examples of FGMs ([4,5] and references quoted therein). Although several 
aspects of FGMs have been reviewed comprehensively ([6-8] and references quoted 
therein) only few investigations have been made to assess the role of the dislocations 



With the exception of [1], the classical theory of linear elasticity was routinely 
utilized to calculate the elastic fields produced by defects (dislocations and disclina- 
tions) in FGMs. However, classical continuum theories are scale invariant in which 
no intrinsic length appears and so fails when one attempts to explain the nano-scale 
phenomena near defects. As a result, elastic singularities are present in these solu- 
tions and size effects which dominate in small volumes cannot be captured. These 
undesirable features are removed within the second strain gradient elasticity formu- 
lation presented in this paper. The solutions obtained herein by using the Fourier 
transform technique and a displacement formulation are reduced to the correspond- 
ing expressions of classical elasticity and first gradient elasticity as, for example, 
were obtained in [1] through the use of the stress function approach. As a result, 
analytical expressions for the displacement field (in addition to those for the stresses 
and strains) are derived. The extra dividend is the derivation of non-singular ex- 
pressions for the double stresses which diverge near the dislocation core in the first 
order strain gradient theory. 



We consider a screw dislocation with Burgers vector b = (0,0,62) in an infinite 
medium with a varying shear modulus = /x(y) = /ioe^'*^ (a ^ 0) in the frame- 
work of classical elasticity. This is a problem of anti-plane shear with the only 
non- vanishing component of displacement u^{x, y) satisfying the displacement equi- 
librium equation 



where denotes the Laplacian. Using the substitution v^^ = e we obtain 



in FGMs. 



2 Classical Solution 




(1) 




(2) 



which by means of the Fourier transform 




where i 




is reduced to the following ordinary differential equation 




(3) 



Since is finite everywhere, we arrive at 




{y > 0), 
(y < 0), 
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where the two unknown functions, A{s) and B{s), are constants with respect to y. 
In view of the present dislocation configuration, we have 



u°(2:,0+) -n°(x,0~) = (5) 



where e^y is the classical strain, and H{—x) is the Heaviside step function. Taking 
the Fourier transform of the above conditions, we have 

4ix,0+) - n°(x,0-) = h [TTdis) + -J ; ^(s,0+) = ^(s,0-), (6) 

and, thus, the unknown functions A{s) and B{s) are determined as 

A{s) = + — , B{s) = - — -Trdjs). (7) 



It follows that 



r,0 



6, e-^y 



2sVs2 + a2 2s 
and by taking the inverse Fourier transform, i.e. 



+ sgn(y) — e-lfl^^ - vr 5{s) H{-y) e^^ 



(8) 



oo 



we obtain (in view of the symmetry of the integral) the final expression 

n° = ^ e-^y r ^^k± e-\y\^ds (9) 
271" Jo sVs^ + a^ 

- sgn(,) |l e-- /"^ ^ e-l^lv^ds - ^ 
27r Jo s 2 

Next we note that in small strain theory the (compatible) total strain efj may be 
written as 

^Ij = {'^/2){ui,j + Uj^i) = Eij + efj, 

where e^j and ef^ denote the usual (incompatible) elastic and plastic strains, respec- 
tively. It follows that 

- sgn{y) ^ e-^y / cos (sx) e-l^^l^^^ds, 
47r Jo 

^OT ^ 1 M = ^ -ay /""^ ^sin(sx) 

2 5y 47r io Vi^T^ 

^ /"°° sin(sx) n 

With the help of the identities 

r '-^ds = ^sgn(x), rJ^e-\y\^ds = ^ar), 
Jo s 2 Jo Vs^ + a^ 



where r = x"^ + y'^ and denotes the modified Bessel function of the second 
kind and of order n, the integrals appearing in Eq. (10) can readily be evaluated to 
give 



-OT 



47r 



-ay 



aK^iar) Ki(ar) 

r 



'zy 



iTT r 2 

(11) 



The last term in Eq. (11)2 which is singular on the half-plane y = and x < 0, 
corresponds to the plastic strain e^y = bz6{y)H{—x)/2. The other term on the 
right hand side of Eq. (11)2 may thus be regarded as the elastic strain. Using the 
constitutive law, a^- = (i = x,y), the stresses read 



27r 



ay 



aKo{ar) Ki{ar) 



a 



2-K 



,ay 



ax 



— Kiiar). 



(12) 



which are the same as those earlier obtained in [1] by the stress function approach, 
and which are singular at the dislocation line. 



3 Strain gradient elasticity solution 

Within a simplified theory of linearized anisotropic theory of second strain gradient 
elasticity proposed in [9] (for a corresponding form of first strain gradient elasticity 
and a robust method for solutions of corresponding boundary value problems, the 
reader may consult [10,11]), the strain energy density has the form 

— 2 ^ijkl^ij^kl ~^ 2 ^ Cij-mn^mn,k^ij,k ~^ 2 ^ Cijfnn^mn,kl^ij,kh (1^) 

where Sij is the elastic strain tensor, i and i' are internal lengths, and Cijki is the 
stiffness tensor of the form 

Cijki = \{x)5ij6ki + /i(a;) {6ik5ji + SjAi) , 

with the Lame constants \{x) and /i(£c) being given functions of the spatial coordi- 
nates. The corresponding expressions for the elastic- like first order stress (f^) and 
the higher-order double (rj^fe) and triple (Tij^i) stresses are given by 



^ij ■ Q — CijlslEklj Tij}^ . — — ^ Cijmn^mn,kj Tijkl ■ — o — ^ Cijmn^mn,kl ■ 

OEij OEij^k C'£ij,kl 



while the Cauchy stress aij (note that in the notation of [1] this was denoted by a^j 
and termed total stress) satisfies, in the absence of body forces, the usual equilibrium 
equation 

'-^ijj ~ '^ijk,kj ~l~ 'Tijkl,klj ~ 0- (1^) 

For the present case of anti-plane shear we have 
afj = 2 fie^j, Tzjk = 2t^ tJ'ezj,k, t^jm = 2 i'^ fJ.ezj,ki; (i, k,l = x,y). 

For an exponentially graded material in the y-direction, i.e. n = fi{y) = /ioe^"^^, it 
follows from the above relations and definitions that the governing equation for the 
gradient displacement Uz (recall the Ru-Aifantis theorem [10,11]) reads 



l-c?(v^ + 2a| 



l-C2(v2 + 2a^ 



Uz = ul, (15) 
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where c\ + = cfc^ = and ti^ denotes the classical elasticity solution 
discussed in the previous section [Eq. (9)]. 

As before, by the substitution Uz = w e""^, we obtain 



[l-cl (V^-a^)] [l-cl (V^-a^)] 



w = w 



(16) 



where is given in Section 2. Use of the two dimensional Fourier transform yields 
the algebraic equation 



[l + cl [s'^ + f + a')] [l + cl{s^ + t' + a')] 



w = w 



(17) 



where 



w = d{d{w; X ^ s};y ^ t} 



oo roo 
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and 
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Then the inverse Fourier transform gives 
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a + it 



a + it 



ia + t cf — C2 t^ + k\ cf — cl t"^ + kI 



(18) 

I e''^''-'+^yUsdt 

I e'^'^'+^yUsdi 
e'^ydL 



where = ya^ + l/c^. If we integrate out the variable s, and use the integral 
relations 



1 



2vr s(s2 + fc2) 



c^s = ^ (l -e I'^^l) sgn(x), 



cos(ty)^^^^^_|,^l 



/o t2 + ifc2 2\k\ 

as well as the symmetry properties of these integrals, can finally be expressed in 
terms of sine and cosine integrals as follows 
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where is the classical solution given by Eq. (9). 

It is easily seen that this expression for a — t- coincides with earlier results 
obtained by the third author and co-workers for homogeneous media and amended 
in [12]. It also turns out that for i = and £' = the result reduces to Eq.(9). 
When X —7- 0, the displacement field can be expressed in an explicit form 



.2 _ „2 
-1 '^2 



cl sgn(y) e '"^'^l + cl sgn(y) e 



-i^2\y\ 



-K2|S/| 



(20) 



It is worth noting that the classical displacement n^(0, y) has an abrupt jump at the 
dislocation line y = 0, while the gradient solution of Eq. (20), is smooth there. For 
a fixed value of o, the larger the ratio C2/C1 is, the smoother the solution becomes. 
It also turns out that the strains are given by 
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3gn(2;) e-l^l\/*^ + 2/7(-x)l dt. (22) 



where e^x ^zy the classical (elastic) strains given in Section 2. It is seen 
that ejj, does not contain a plastic part, while e^^^ is decomposed into the elastic 



and plastic strains, i.e. 
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0, and 
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Air ci 
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- [-clKiKiimr) + clK2Ki{K2r)] , ef^ = ej^ - F^y 



It is also follows that the expression for the plastic strain e^(0,y) is given by the 
simple formula 
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The lower-order elastic-like stresses and higher-order double stresses are given 
by the expressions: 
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where a^^ and a^y are the classical stresses given by Eqs. (12), and 
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It is seen from the above expressions that af^ is still symmetric with respect to the 
plane x = 0, while afy has lost symmetry with respect to plane y = 0. Moreover, 
in contrast to homogeneous medium, T(^zx)x / ~T{zy)y More details and for this 
problem their physical implications to possible improvements of designing FGMs 
and the expressions for the triple stresses will be given in a forthcoming publication. 
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